The next to leading order effective potential in the 2+1 dimensional 
Nambu-Jona-Lasinio model at finite temperature. 
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The finite temperature effective potential in the 2+1 dimensional Nambu-Jona-Lasinio model 
is constructed up to the next to leading order in the large N expansion, where N is the number 
of flavors in the model. The distinctive feature of the analysis is an inclusion of an additional 
scalar field, which allows us to circumvent the well known, and otherwise unavoidable problem 
with the imaginary contribution to the effective potential. In accordance with the Mermin-Wagner- 
Coleman theorem, applied to the dimensionally reduced subsystem of the zero Matsubara modes of 
the composite boson fields, the finite temperature effective potential reveals a global minimum at 
the zero of the composite order parameter. This allows us to conclude that the continuous global 
symmetry of the NJL model is not broken for any arbitrarily small, finite temperature. 
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I. INTRODUCTION 



The Nambu-Jona-Lasinio (NJL) model plays a crucial role in our understanding of quantum field theory Despite 
its non-renormalizability in 3+1 dimensions, the NJL model has already found many applications in nuclear and high 
energy physics |,|| . 

£T) , The 2+1 dimensional NJL model has also been studied in the large N limit This model is know to be non- 

£NJ ' renormalizable when expanded in the four-fermi coupling. Nevertheless, it is renormalizable in the 1/N perturbation 
theory PLQ. The effective potential of this model has been studied at zero temperature up to the next to leading 
order in 1/N by carrying out the Hubbard- Stratonovich transformation H,|). When the fluctuations of the composite 
scalars are included the effective potential develops an imaginary part. At finite temperature this complicates the 
analysis of the symmetry breaking and restoration. 
i-SZh ' In this paper, we study the effective potential to the next to leading order in the large N expansion for the three- 
dimensional model at zero and finite temperatures. The problem of obtaining the effective potential is an old one. 
There exist well developed general methods for constructing effective potentials for elementary [§-|Io| as well as for 
D [ composite fields |ll| (see also |l2j where a finite temperature version of the latter is discussed). 

• These methods, however, cannot be applied to some specific models in a straightforward manner. For example, 
| should one try to construct the effective potential for a scalar theory with a double-well tree level potential, he 
would find at the leading order an imaginary contribution |^,|l3| [is]]. The latter is due to a tachionic vacuum at an 
intermediate stage in the calculations. In the scalar 4> 4 theory, a method for dealing with this obstacle was found long 
ago pq-p^|. The way around the problem was the introduction of an additional scalar field. 

While studying the NJL model, one usually introduces some auxiliary scalar fields that could be interpreted as 
composites built from fermions ||,f7| . In terms of these new fields, the problem of constructing the effective potential 
is very much the same as in pure scalar theories except that the propagators of the composite particles are much more 
complicated 

The leading order effective potential in the NJL model, which plays the same role as a tree level potential in scalar 
theories, may also be double-well in shape. In fact, this happens as soon as the four-fermi coupling constant becomes 
larger than some critical value. Then, an attempt to go beyond the leading order inevitably results in an imaginary 
contribution to the effective potential due to the same reasons as in the pure scalar case Q. 

Often there is no reason to doubt the leading order result, so this lack of a straightforward method for calculating 
the next to leading order corrections may seem to be of no interest. There are, however, several instances when 
the knowledge of the next to leading order terms becomes of prime importance. This happens, for example, when 
one studies the finite temperature 2+1 dimensional NJL model which enjoys a continuous global symmetry. The 
leading order calculation of the effective potential indicates that the symmetry remains broken for a range of non-zero 
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temperatures. This result seem to contradict the predictions of the Mermin-Wagner-Coleman (MWC) theorem p0[ . 
The argument is well known. The infrared region of the system is dominated by the zero Matsubara modes of the 
boson fields. The latter effectively live in a 1+1 dimensional space. Thus, by appealing to the above mentioned 
theorem, one concludes that the symmetry breaking is forbidden |pi|fffi| . 

We find it surprising that there are no explicit calculations of the next to leading order effective potential at finite 
temperature in the NJL model in 2+1 dimensions to demonstrate that the symmetry is restored for an arbitrarily small 
T =/= 0. In this paper we perform such calculations and show that the global minimum of the potential does indeed 
occur at the origin or, in other words, that the symmetry is not broken at any finite temperature. The fluctuations 
of the composite scalar particle, which is the NG boson at T — 0, play a crucial role in the symmetry restoration. 

In some recent studies the finite temperature condensate f2^,|2J] and the equation of state of NJL models have 
been computed to the next to leading order in 1/N. It is worthwhile to note that the problem of the imaginary 
contributions was appreciated in pq] . There it was noticed that the effective potential of the NJL model (with a 
non-zero bare mass for fcrmions) is real at the minimum field configuration. This is similar to the result obtained by 
Dolan and Jackiw for theories with fundamental scalar fields |l|. Even though the effective potential is real at the 
minimum, the problem still exists in general. In this paper we suggest a method that resolves the imaginary part 
problem in the NJL model in 2+1 dimensions at finite temperature (and even outside the minimum of the potential). 
Its generalization may also work in the 3+1 dimensional case at finite temperature. 

Finally, it seems that there is no clear understanding of what happens in the 2+1 dimensional NJL model as the 
temperature approaches zero. Since the MWC theorem is believed to work for arbitrarily small temperatures, there 
is no symmetry breaking for any non-zero T. At T=0, the continuous symmetry is broken. Then, a natural question 
arises. How does the effective potential in the model behave as the temperature approaches zero? We answer this 
question below. 



II. ZERO TEMPERATURE CASE 



Let us start with the analysis of the NJL model in 3-dimcnsional Euclidean space at zero temperature. To apply 
the large N expansion scheme, we consider the model containing N identical flavors. The Lagrangian density reads 
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J NJL 
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(1) 



where we have chosen four-component spinors and the 4x4 Dirac 7-matrices are all antihermitian. By introducing 
("integrating in") the auxiliary boson fields a = (G/N)^^ and it = (G/N)^>ij5^> , we obtain the equivalent model 
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The mentioned equivalence comes from the Euler-Lagrange equations of motion for the boson fields. As is easy 
to check, the latter Lagrangian density enjoys the same continuous chiral symmetry as that in Eq. (Q). The only 
difference is that the integrated in boson fields are also involved in symmetry transformations. 
After integrating out the fermions, we are left with the following effective action for the bosons: 
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This action may be expanded in inverse powers of the cutoff A in the following manner, 

~6 ' 



S 



(o) 

seal 



N 




A 






/ (fx 




H) 






-2 





J , P 3 , P 4 



8A 



O 



(3) 



(4) 



where p = v a 2 + n 2 is chiral invariant, and the dimensionless coupling constant is defined as g = GA/tt. In the kinetic 
part of the action, we keep only the terms that contain up to two derivatives, i. e., terms quadratic in momenta, 
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The straightforward calculation of the infrared asymptotic behavior of the boson propagators reveals the following 
expressions for the kinetic coefficients p(|p7[ 
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(6) 



Instead of using these infrared asymptotes, we replace them with constant coefficients and argue that this would give 
a reasonable approximation in the whole range of momenta. The reasoning goes as follows. The composite particles 
in the far infrared region (large distances) should presumably behave as almost point-like objects. Their kinetic 
coefficients would be given by the expressions as in Eq. (^) with a specific value of p = p. As for the ultraviolet region, 
it is known that the leading asymptotic behavior does not depend on p. Therefore, we take a^'s independent of p as 
a much better approximation than that in Eq. (^). Later, when we match the zero and the finite temperature results, 
we will see that there is another advantage to the above approximation. 

Following the standard analysis applied for scalar theories [|l6|,|l8| , we introduce an additional auxiliary scalar field 
X into our bosonic action, 
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where p = A(g — l)/g is the characteristic scale of chiral symmetry breaking. Below we assume that p <C A. As 
we mentioned in the Introduction, the model at hand is renormalizable in the 1/N expansion. In spite of that, it is 
more convenient for us to use the bare quantities throughout the paper without performing renormalizations. That 
is sufficient for our purposes since we are only interested in analyzing the long distance behavior of the theory, which 
is captured quite well by a model with a large but finite physical cutoff. 

In calculations, we have to keep all contributions up to the order 1/A. However, it will be more convenient to keep 
exact results as they appear after integrations, without continually performing the expansions. Obviously, this will 
not affect any of our results since all the interesting features of the effective potential, as we shall see, occur at values 
of a much smaller than the cutoff A. 

To construct the effective potential, we apply the standard technique jlTJ to the action in (Q). After performing 
rather simple calculations up to the next to leading order of the large N expansion, we obtain the following effective 
potential: 
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These parameters define the bosonic fluctuations. Since their propagators do not have the canonical form, and 

r(0) 



Mff should not be confused with the sigma and pion masses 
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Finally, to get the effective action for the a field, we have to eliminate the x ne ld by the saddle point equation 
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We solve this equation numerically to obtain x = x( a )- Then, we substitute it into Eq. (@) to get the effective 
potential as a function of a. The minimum of the obtained effective potential satisfies the following gap equation, 
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Prior to presenting the numerical results, let us analyze the saddle-point equation for the x held given in Eq. (|1] 
It is easy to see that the left hand side of this equation is a monotonically increasing function of X- We also observe 
that this function takes real values only for x > ~° "■> an d so, it is bounded below with the minimum at x = ~ a ■ If 
this minimum value is positive, the left hand side remains positive in the whole range x ~ a i an( I f ne saddle point 
equation does not have a real solution. Therefore, the existence of a real solution x = x( a ) requires the following 
condition to be satisfied, 
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This condition may be rewritten in a much more transparent way. Indeed, if we choose a large enough value of N 
the condition in Eq. (|l^) does not hold for a = 0. By combining this with the fact that the function in the left hand 
side of Eq. (^) decreases with increasing a (this is true everywhere except for a small unimportant region of order 
/j, 3 /(NA) 2 around a — 0), we conclude that the above condition is equivalent to a > a cr , where 
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For any value of a greater than a cr , we can numerically solve the saddle point equation (|l l|) . Then substituting 
the solution \ — x(cr) into Eq. (g), we finally obtain the effective potential as a function of a alone. Our numerical 
analysis shows that the potential is a monotonically increasing function in the region a > a cr for a very wide range of 
input parameters, which are the cutoff (A), the temperature (T), the number of flavors (N), the kinetic coefficients 
{a£\ and a^), and finally the scale fi = A(g — l)/<?, which we put to 1 without loss of generality. In Figure 1, we 
present the plot of the effective potential in the region a > a cr obtained for a specific set of the parameters. 

Obviously, the straightforward extension of our calculations to the region < a < a cr would result in a non-zero 
imaginary part in the effective potential. Such a situation looks rather disappointing. It has been claimed, though, 
that the Maxwell construction has to be used for the region < a < a cr [28[. The latter solution of the problem 
presumably corresponds to the flat potential in the "forbidden" range of a as is plotted in Figure 1. We would like 
also to mention Ref. p^j where a strong argument in support of the flat potential was given. This behavior of the 
potential is interpreted as a signal for the phase with the broken symmetry. 



III. NON ZERO TEMPERATURE 



Now, let us study the NJL model at finite temperature. The partition function is given by the following path 
integral 



Z = j ( J[ d^dV J 



e - bNJL , (15) 



where the functional measure is defined in a way consistent with the standard anti-periodic boundary conditions for 
the fermion field. The Euclidean "action" that appears in the partition function is given by 
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where (3 = l/T is the inverse temperature. 

As in the case of zero temperature, we introduce two auxiliary scalar fields, so that the partition function is given 
by a path integral over both the fermion and the boson fields. The corresponding action reads 
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We are interested in the effective action for the cr-field. In order to obtain it, first, we calculate the expression for 
the partition function in terms of the scalar fields alone by performing the Gaussian integral over the fermions. As a 
result, we are left with the following "exact" action for the scalars, 
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and the partition function is given by the integral over the scalar fields, 
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Following closely our zero temperature calculation, we approximate the action in Eq. ( |l8| ) by a truncated derivative 
expansion (which in momentum space corresponds to the expansion in powers of momenta up to the second order). 
The approximate action is 
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where we use the same notation as in the previous section. The kinetic term looks somewhat different from that at 
zero temperature, 
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where the coefficient a's and o's are now not equal. This is due to the heat bath which always breaks Lorentz 
invariance. In particular, the infrared asymptotic behavior of the kinetic term leads to the following coefficients, 
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As in the zero temperature case, we do not use these expressions later in the calculation, but switch to the same 
approximation as before. Namely, we assume that the a's and b's are independent of p. Later we give arguments to 
substantiate this claim. 

After "integrating in" the auxiliary x field, we arrive at the following action, 
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This is used for the construction of the effective potential. As in the case of zero temperature, we obtain the result of 
interest in the next to leading order of the large N expansion: 
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After carrying out the summation over Matsubara frequencies and the integration over momenta (with the sharp 
cutoff at A), we arrive at the result 
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Despite the fact that the calculations were simplified due to the nontrivial assumption on the form of the kinetic part 
of the boson action, we argue that the result should be reliable. Indeed, the essential details of both the infrared and 
the ultraviolet regions seem to be captured if the kinetic coefficients are chosen to be independent of a. 

The only property that we require of the boson propagators in the ultraviolet region is a soft dependence on the 
field <7. On the one hand, this resembles the known asymptotic behavior of the composite boson propagators |fl^9[, 
and, on the other, this allows us to avoid the problem of an unphysical sensitivity on the details of the ultraviolet 
region. 

As for the infrared behavior of the propagators, the argument is somewhat trickier. We have to distinguish between 
small and large values of the field a. Obviously, the most interesting region lies close to the origin, since it is there 
where one would expect the essential changes in the effective potential when the chiral symmetry is restored. So, 
suppose that we consider a < T. Then for \p\ < T, the kinetic coefficients approach constants of order 1/T. Since 
only the zero Matsubara mode is relevant in this case, we can also ignore the fact that b a ,n — > 0. On the other hand, 
if we consider the region a > T, then, strictly speaking, we have to take into account any dependence of the a's and 
b's on the value of the field a. Fortunately, there is no reason to believe that a mild dependence of the propagators 
on a in this region could change any of our qualitative results. 

Now, let us mention that the naive limit T — > of our finite temperature potential given in Eq. ( ^0| ) coincides with 
the effective potential obtained in the previous section for T = if we assume that a a ,b a — > a§ and 0^,6^ — ► ai°\ 

To eliminate the auxiliary field Xj we have to solve the saddle point equation, as was done in the case of zero 
temperature. This saddle point equation is obtained by taking the derivative of the effective potential (|3^) with 
respect to x : 
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Let us analyze this equation. As in the case of T = 0, the left hand side of Eq. ( pl| ) is a monotonically increasing 
function of \- Unlike the former, though, the left hand side of eq. ( |3l| ) is not bounded below, but approaches negative 
infinity as x ~~ * ~ a ~ 2Tln[l + cxp(— er/T)], so that M w — > 0. This simple property has far reaching consequences. 
The most important among them is that the saddle point equation has a real solution for all values of a ! This is quite 
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remarkable if we notice that Eq. fl31| ) in the "naive" limit T — > coincides with its counterpart that was obtained for 
T = 0. 

The saddle point equation for T — 0, as we saw, does not have a solution for < a < <J cr , while the finite 
temperature equation allows a real solution for all values of a. How could this happen? The answer is quite simple, 
and involves the way the limit T — > is taken. In the "naive" limit, we assume that M n and M a become much larger 
than the value of T as we approach zero temperature. On the other hand, at any small but finite temperature, we find 
that the solution to the saddle point equation requires that M w ~ T 2 b-x exp(— NAf(a)/T), where f(a) is a positive 
function of a for < a < cr cr (T). 

Therefore, we observe that the T = quantum field theory corresponds to the "naive" limit when T — > is taken 
before solving the saddle point equation. We get a very different result by taking the limit in the other order, i.e., 
solving the saddle point equation first and then taking the value of T to zero. 

The finite temperature gap equation is given by 




= 0, 



(32) 



and admits the trivial solution a = 0. For T = 0, this was not the case. Actually, the existence of the trivial solution 
is a necessary condition for concluding that the symmetry is restored, but not a sufficient condition. Below, while 
presenting our numerical results, we will see that this a — solution is the only solution to the gap equation, and 
that it corresponds to the global minimum of the effective potential. 

Our numerical analysis begins with a solution of the saddle point equation (|3l|). While solving it, we notice that 
there are two qualitatively different regions of the field a. For small enough values of the field a < a cr (T), where 
the T = equation would not have a solution, we can give a very good analytical estimate for the finite temperature 
solution, 



where 
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and where f(a) is a positive function of a for < a < a cr (T). After a reaches the value a cr (T), the solution x(c) 
sharply changes its behavior and starts to follow closely the solution found for T = 0. The results of our numerical 
solution to the saddle point equation for four different temperatures are plotted in Figure 2. 

By substituting the solution \ — x(°") of the saddle point equation into Eq. (fujl), we obtain the effective potential 
as a function of a. The corresponding plots for the four different temperatures are given in Figure 3. As we see, the 
effective potential reveals the behavior which was expected: its global minimum lies at the origin. The latter means 
that the chiral symmetry is not broken. 

By comparing the zero and finite but small temperature effective potentials, we find that they follow very closely 
each other in the range of large enough values of the field a > a cr . 

Another property of the obtained effective potential which is worth noting is its convexity. This is a feature that 
one expects for the effective potential based on general arguments ||l]] . To see that our effective potential is indeed 
convex, we have to show that its derivative is non-negative in the whole range of field a. For this purpose, in Figure 4, 
we plot the left hand side of the gap equation (|3^). Up to an unessential factor, it is equal to the derivative of the 
potential as a function of a. As we see, the latter is everywhere non-negative. 

Before concluding this section, we have to note that all the qualitative results, such as the existence of one global 
minimum at a = 0, the convex character of the potential, and the two-branched character of the solution to the saddle 
point equation, persist for a rather wide choice of the input parameters. 



IV. CONCLUSION 



In this paper we have studied the effective potential in the three-dimensional NJL model at finite and zero tem- 
peratures up to the next to leading order in the large N expansion. The distinctive feature of our analysis is the 
introduction of an additional scalar field. This allows us to circumvent the well known, and otherwise unavoidable, 
problem of the imaginary contributions to the finite temperature effective potential. 
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In agreement with the MWC theorem, the analysis of the effective potential in the three-dimensional NJL model at 
finite temperature does not reveal breaking of the continuous chiral symmetry. This remains true for an arbitrarily 
small value of the temperature. The leading order effective potential that indicates symmetry breaking changes 
drastically in the next to leading order. The essential role is played by the pion-like particle, which would have 
become the NG boson were the symmetry broken. 

We present the numerical results that clarify and confirm our analytical findings in Figures 1-4. The zero temper- 
ature effective potential is given in Figure 1. Since it is not defined in the range of fields < a < er cr , we extended it 
with the Maxwell construction to the forbidden region. More rigorous arguments in support of such an extension can 
be found in p9fl . Figures 2 through 4 present the results at finite temperature. 

By comparing our analytical and numerical results for both T = and T ^ 0, we come to the conclusion that the 
zero temperature field theory can be defined in two different ways. One definition comes from the standard T = 
quantum field theory. The other is obtained by the limiting procedure T — ► from the Matsubara imaginary time 
formalism. 
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